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Abstract 



^-H Secure communications can be impeded by eavesdroppers in conventional relay systems. This paper 

o 

^NJ proposes cooperative jamming strategies for two-hop relay networks where the eavesdropper can wiretap 



the relay channels in both hops. In these approaches, the normally inactive nodes in the relay network 
can be used as cooperative jamming sources to confuse the eavesdropper Linear precoding schemes 
are investigated for two scenarios where single or multiple data streams are transmitted via a decode- 
and-forward (DF) relay, under the assumption that global channel state information (CSI) is available. 



jy! For the case of single data stream transmission, we derive closed-form jamming beamformers and 

I the corresponding optimal power allocation. Generalized singular value decomposition (GSVD)-based 

,__( secure relaying schemes are proposed for the transmission of multiple data streams. The optimal power 

allocation is found for the GSVD relaying scheme via geometric programming. Based on this result, a 

GSVD-based cooperative jamming scheme is proposed that shows significant improvement in terms of 

secrecy rate compared to the approach without jamming. Furthermore, the case involving an eavesdropper 

with unknown CSI is also investigated in this paper. Simulation results show that the secrecy rate is 

^~~^ dramatically increased when inactive nodes in the relay network participate in cooperative jamming. 

> 

; ^ Index Terms 

X 

C^ Wiretap channel, secrecy, relay networks, physical layer security, jamming, interference. 
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I. Introduction 

Security is an important concern in wireless networks due to their vulnerability to eavesdropping. 
Traditionally, security is viewed as an issue addressed above the physical (PHY) layer, and all widely 
used cryptographic protocols are designed and implemented assuming the physical layer has already been 
established and provides an error-free link (V\. However, higher-layer key distribution and management 
may be difficult to implement and vulnerable to attack in complex environments such as ad-hoc or relay 
networks, in which transceivers may join or leave randomly Q, Q. Therefore, there has recently been 
considerable interest in physical layer security, which explores the characteristics of the wireless channel 
to improve wireless transmission security. 

The theoretical basis of this area was laid by Wyner, who introduced the wiretap channel and demon- 
strated that when the eavesdropper's channel is a degraded version of the channel of the legitimate 
receiver, the transmitter can send secret messages to the destination while keeping the eavesdropper from 
learning anything about the message |j4|. The notion of secrecy capacity was introduced and defined as 
the maximum achievable transmission rate of confidential information from the source to its intended 
receiver. Later, Csiszar and Komer generalized Wyner's approach by considering the transmission of secret 
messages over broadcast channels [5]. Recently, considerable research has examined secrecy in wiretap 
channels with multiple antennas |[6|-|[T4|. In particular, the secrecy capacity of the multiple-input multiple- 



output (MIMO) wiretap channel has been fully characterized in |10|, |11|. With the additional degrees 



of freedom provided by multi-antenna systems, transmitters can generate artificial noise to degrade the 



channel condition of the eavesdropper while maintaining little interference to legitimate users |13|-|16|. 



As a natural extension, approaches for physical layer security have also been investigated in cooperative 



relaying networks 1 17|-p2|. In these cases, relays or even destinations can be used as helpers to provide 



jamming signals to confuse the eavesdropper. This approach is often referred to as cooperative jamming. 



In 1 20 1, a noise-forwarding strategy is introduced for a four-terminal relay-eavesdropper channel where 
the full-duplex relay sends codewords independent of the secret message to confuse the eavesdropper. 
A two-stage cooperative jamming protocol is investigated in |14], where multiple relay nodes act as an 
extension of the single-antenna source node. In this work, the "relays" only play the role of a helper and 
do not relay the information signals. In [211 , three cooperative schemes are proposed for a single-antenna 
relay network, and the corresponding relay weights and power allocation strategy are derived to enhance 
the secrecy for the second hop. An optimal beamforming design for decode-and-forward (DF) relays is 



investigated in |22|, but only the scenario where the eavesdropper wiretaps just the link between the relay 



and destination is considered. 

Unlike the aforementioned work, this paper proposes cooperative jamming strategies for a half-duplex 
two-hop wireless MIMO relay system in which the eavesdropper can wiretap the channels during both 
transmission phases. Cases involving both single and multiple data stream transmissions are investigated. 
Due to the lack of "outer" helpers, the source, relay and destination must rely on themselves for jamming 
support. This approach guarantees that the eavesdropper is jammed whether it is close to the source or 
the destination. In the proposed cooperative jamming strategies, the source and the destination nodes 
act as temporary helpers to transmit jamming signals during the transmission phases in which they 
are normally inactive. We define two types of cooperative jamming schemes, full cooperative jamming 
(FCJ) and partial cooperative jamming (PCJ), depending on whether or not both the transmitter and the 
temporary helper transmit jamming signals at the same time. 

We focus on the design of linear precoding schemes throughout the paper, and begin with a simple 
scenario where the relay has only a single antenna. In this case, we investigate the joint design of the 
jamming beamformer and the power allocation for two optimization problems: (1) maximizing the secrecy 
rate with certain power constraints, and (2) minimizing the transmit power with a fixed target secrecy 
rate. Since a joint optimization of the beamformers and power allocation is in general intractable even 
if global CSI is available, we use a suboptimal zero-forcing constraint that the jamming and information 
signals lie in orthogonal subspaces when received by the legitimate nodes, and we derive closed-form 
expressions for the jamming beamformers. Based on these results, we find the optimal solution for the 
power allocation by utilizing the method of geometric programming (GP). Then we expand the scope 
to study the scenario where all nodes have multiple antennas, and multiple data streams are transmitted 
via the relay. A generalized singular value decomposition (GSVD)-based cooperative jamming scheme 
is proposed and the corresponding power allocation strategy is discussed. Unlike the single data stream 
case that uses a zero-forcing constraint, the cooperative GSVD-based jamming method will not in general 
produce jamming signals that are orthogonal to the desired signal. 

Another important consideration is the availability of the eavesdropper's CSI. If the CSI of the 
eavesdropper is known, (for example, if the eavesdropper is another active user in the wireless network), 
the transmitter can optimize its beamformer to enhance the information transmission to intended nodes 
while suppressing or even eliminating the leakage to eavesdroppers. However, in some cases {e.g., passive 
eavesdroppers), it is impractical to assume known CSI for the eavesdroppers. Since the secrecy rate can 
not be optimized without knowledge of the eavesdropper's CSI, we will follow the approach of |fT5l, 



1 25 1, where the transmitter first allocates part of its resources to guarantee a fixed target rate, and 



then uses the remaining resources to jam the eavesdropper. 

The organization of the paper is as follows. Section[II]describes the system model considered throughout 



the paper. In Section III the cooperative jamming schemes, including the jamming beamformer design 
and power allocation, is investigated when the eavesdropper's CSI is known. Both single and multiple data 
stream transmissions are considered in this section. Secure relaying under the assumption of unknown 



eavesdropper's CSI is studied in Section IV The performance of the proposed cooperative jamming 
schemes are discussed in Section |Vj and conclusions are drawn in Section |Vl] 

The following notation is used in the paper: IE{-} denotes expectation, (•)^ the matrix transpose and 
(•)^ the Hermitian transpose. || • || represents the Euclidean norm, | • | is the absolute value, [x]^ denotes 
max{x,0}, tr(-) is the trace operator, J\f{-) represents the null space, and I is an identity matrix of 
appropriate dimension. 



II. System Model 

We consider a two-phase four-terminal relay system composed of a source (Alice), a destination (Bob), 
a DF relay node and an eavesdropper (Eve), as shown in Fig. [T] The message from Alice is uniformly 
distributed over the message set W = {1, • • • , 2"^}, where R denotes the source rate in bits per channel 
use. The confidential message is randomly mapped to a length-n source codeword z" € Z^ and the Relay 
encoder maps its received signal to codeword z" G Z^, where ^^ and Z^ are length-n input alphabets. 
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Fig. 1. Relay scenario. 



All nodes are assumed to be half-duplex, i.e. a two-hop time division multiple access system is 
considered. Alice transmits in the first phase while the relay listens, and relay transmits in the second 
phase. We assume there is no direct communication link between Alice and Bob, except perhaps for 



some low-rate control or channel state information, and thus Alice and Bob must rely on two-phase 
transmissions through the relay. This is a reasonable assumption in the type of scenarios where relaying 
is used, where direct high-rate communication is too "expensive" in terms of the given power constraints, 
but low-rate control information can still be exchanged [19]. When Alice transmits a jamming signal, 
however, its impact on Bob's received signal must be taken into account. All nodes in general have 
multiple antennas. The number of antennas possessed by Alice, Bob, the Relay and Eve are denoted 
by Na, Nh, Nr and N^., respectively. In part of the paper, we will explicitly consider scenarios where 
the Relay has only a single antenna. We restrict attention to scenarios where all nodes (including the 
eavesdropper) employ linear preceding and receive beamforming. 

A. Relay Transmission 

In the first phase, Alice transmits the information signal to the Relay. Both the Relay and Eve will 
receive the signal as 

Yr = HarTaZa + n^ (1) 

Yel = UaeTaZa + Uel (2) 

where Zq is the information signal vector transmitted by Alice, T^ G c^"^*^ {1 < k < s) is the transmit 
beamformer used by Alice, and we assume m = ranklHar}, n = ran/c{H,.b}, s = min(?n,, n) and 
k represents the number of data streams to be transmitted. The terms n^ and ngi represent naturally 
occurring noise at the Relay and Eve, respectively. For simplicity, we assume that the noise vectors at all 
nodes are Gaussian with covariance cr^I. In general, Hjj (hjj) represents the channel matrix from node 
i to j, with i,j £ {a, 6, e,r} denoting which of the four terminals is involved. These channel matrices 
are fixed over both hops. The signal received by Bob and Eve in the second transmission phase can be 
expressed as 

Yb = HrftT,.z,. + rif, (3) 

ye2 = HreTrZr + Ta.e2 (4) 

where Zj. is the signal vector transmitted by the Relay, T^ G C^'' ^ '^ is the transmit beamformer used by 
the Relay, and n;,, ne2 represent the noise vectors at Bob and Eve. There is a transmit power constraint P 
on both phases, i.e., E{z^Za} < P and Efz^z^} < P. We assume a repetition-coding scheme, where z,, 
is simply a scaled version of z^. In particular, we assume z^ = D^z and z^ = D^z, where E{zz^} = I 
and Da, D^ ai^e diagonal power loading matrices that ensure the power constraints are met. 



B. Cooperative Jamming 

In the most general case, the signals transmitted by Alice in the first phase may contain both information 
and jamming signals, and Bob may also transmit jamming signals at the same time. Thus the signals 
received by the Relay and Eve in the first phase will be given by 

y, = Har(TaZa + T'^z'J + Hfc^T'fez'f, + n, (5) 

yel = H„e(TaZ« + TM + HbeT'feZ; + Ilel (6) 

where z^ and z'^ are jamming signal vectors transmitted by Alice and Bob, respectively, and T'^ and 
T'f, are the corresponding transmit beamformers. In this paper, T'^ and T'^ could be chosen to project 
the jamming signals on the subspace orthogonal to the information signals, or they could allow a small 
amount of interference leakage to the legitimate receiver while producing more interference power at 
Eve, as will be discussed when the GSVD-based transmission strategy is used. We refer to the case 
where both t!^ ^ and 2,'^ 7^ as full cooperative jamming (FCJ). If either of them is zero, we refer to 
it as partial cooperative jamming (PCJ). FCJ will not be considered in the scenario where Eve's CSI is 
known, since in this case splitting the power between data and jamming signals at Alice is known to be 
suboptimal. However, when Eve's CSI is known, we will still study the PCJ scheme where Bob uses part 
of the global transmit power to produce jamming signals. When Eve's CSI is not available, FCJ should 



be used, as will be discussed in Section IV 



In phase 2, the signals received by Bob and Eve are given by 

Yb = H,6(T,Z, + T;z;) + nabT'a2<2 + "b (7) 

ye2 = H^e(T,Z, + T'X) + ^aeT'a2Za2 + ^e2 (8) 

where Zj. is the information signal vector of the Relay with transmit beamformer T,., z^ and z'^2 ^^ 
jamming signal vectors transmitted by the Relay and Alice, respectively, and T',. and T^2 ^^ their 
corresponding transmit beamformers. Note that, although there is no direct link for the information 
signal. Bob still sees the jamming signal from Alice. For a global power constraint, we have 

E{zf z, + zf z', + zf z'J < P E{zf z, + zf z; + z'.f z'.J < P. 

We will also investigate scenarios with individual power constraints, i.e. E{z^Za + tI^t!^ < Pa, 
E{zf z'J < Pb, E{z^z, + z;^z;} < P,, and Ejz'^f z^^s} < P„. 



C. Performance Metric 
MIMO wiretap channels have been extensively analyzed in recent work, and the achievable secrecy 



rate has been shown to be 1 10|, |11 1 



Rs = max [Id - /e 



(9) 



where Id is the mutual information from the source to the destination, /e is the mutual information from 
the source to the eavesdropper, and the maximum is taken over all possible input covariance matrices. 
For the half-duplex two-hop relay channel, the achievable secrecy rate was found in |26] to satisfy the 
same expression as in (|9]), where amplify-and-forward, decode-and-forward, and compress-and-forward 
relaying modes were all investigated. Eq. ([9]) was also used as a performance metric to evaluate cooperative 
jamming schemes for half-duplex relay networks in |[2i|. In general, to obtain the maximum secrecy rate, 
one must construct an optimal coding scheme, although potentially suboptimal Gaussian codebooks are 
assumed in ||8|, |21 1, |26|. In Section III we will follow the convention adopted in |21 1, |26| and use Q 
as our metric for evaluating the achievable secrecy rate, assuming Gaussian inputs. Note that Q was 



shown to be valid for both independent and repetition codebooks |26|, although we will only focus on 



repetition coding {e.g. |27|, |28|) at the relay since independent codebooks are expected to result in 
smaller secrecy rates when the encoding schemes and relay protocols are public information p6| . 

The discussion above applies to the cases where the eavesdropper's CSI is known or at least partially 
known {e.g. the case where only statistical channel knowledge is available and ergodic secrecy rate 
is studied ||8|, |29|). However, when the eavesdropper's CSI is completely unavailable, Q may not 
represent an achievable secrecy rate. Some recent progress has been made on finding expressions for the 
achievable secrecy rate in certain scenarios where the eavesdropper's CSI is completely unknown )30| , 
but the derivation of such an expression for the relay network considered here is still an open problem. 
Nonetheless, the difference in the mutual information between the desired receiver and the eavesdropper 
is still a valid metric for evaluating the relative security of competing physical layer approaches. While the 
transmission parameters cannot be chosen to optimize ([9]) when the eavesdropper's CSI is unknown, the 



approach of p5| , p4| ], ||25j] can be followed in which attention is restricted to obtaining a certain desired 
QoS for the legitimate receiver, and then finding a robust strategy for using the remaining resources to 
jam potential eavesdroppers. This is the approach adopted in Section IV with (|9]) as the performance 
metric. 



III. Secure relaying with known ECSI 

In this section, we assume that Eve's CSI (ECSI) is available to the relay network. We will begin 
with the simple case where the Relay is equipped with only a single antenna, then a more complicated 
scenario with a MIMO relay will be investigated. 

A. Single data stream relaying 

We begin by assuming a single-antenna DF relay {N^ = 1), where only one data stream can be 
transmitted via the Relay. Under the PCJ approach, the signals received in each phase can be expressed 
as 

yr = h.artaZa + hferTfeZ'j, + n^ (10) 

yel = HaetaZa + HfeeT'^z'f, + Hel (11) 

and 

yh = hrbZr + ilabT'^z^ + nb (12) 

ye2 = KeZr + H^eT'^z'^ + ne2 (13) 

where E{zf z^} = pa, E{z'/^z'J = pb, E{z^ Zr} = Pr and E{z'/^z'„} = pa2- This is the PCJ form of 
([5])-([8]l with z'^ = and z'^ = 0. Since Nr = 1 in this case, we can design T'^ such that the jamming 
signals aie completely nulled at the Relay, i.e., hbrT'^ = 0. For the transmit beamformer t^ in the first 
phase, we choose the generalized eigenvector of the pencil (1+ 2|h^.har, 1 + %H^Hae) with the largest 
generalized eigenvalue, which achieves the secrecy capacity for the single-hop MISO wiretap channel 



|13|. For the second phase, we design T'„ such that Ha^T'^ is orthogonal to the one-dimensional signal 
subspace spanjhrf,}, so that the jamming does not impact Bob's reception of the information signal. 

1) Maximum secrecy rate with power constraints: Next, we will discuss the design of the jamming 
beamformers and power allocation for maximizing the secrecy rate under both global power constraints 
{Pa + Pfe < ^ in the first phase and pr + Pa2 < P in the second phase) and individual power constraints. 
For a two-hop DF-based relay channel, the mutual information between Alice and Bob through the relay 
link can be written as | |3l| 

1 

-^d = 2 min{log2(l + 7ar),log2(l + 7r6)} (14) 



where | appears because the relay transmission is divided into two stages, and "fij is the SINR at node 
j for the signal from node i. Eve receives data during both phases, and her mutual information is 

1 

/e = - min{log2(l + 7ar),log2(l + 7ae + Ire)}- (15) 

Thus, the secrecy rate can be expressed as 

2 log2 (1+^ +^ ) > lae + Ire < lar < Irb Or -far > max{-frb, lae + Ire] 

foe ,re; ^^-j 

0, otherwise. 

Since the rate of the relay link is limited by the SINR of the inferior phase, for a single data stream 
the transmit power for Alice and the Relay should be adjusted such that ^ar = 7r6 for power efficiency. 
Thus Rs = ^ log2 ,-^^ ^^" N will be used as the objective function in the remainder of this section, as 
a result of the power adjustment. 

We assume Eve uses beamformers Wgi and We2 to receive the signals from Alice and the Relay in 
the first and second phases, respectively: 

W^yel = W^ {YiaetaZa + ^beT^'b'^'b + "el) (17) 

W^ye2 = W^(hreZr + HaeT'^z'^ + ne2), (18) 

and we assume that Eve can compute the beamformers which yield the best SINR, 

Wel = (HfceT'.Q.fc-Tf Hg + a^iy^-tiaeta (19) 

We2 = (H,eT;Q,„,Tf Hf, + a^iy^Ke (20) 

where Q^f,/ = Ejz'^z'j^} and Q^q' = Ejz'jjz'/^}. With the above assumptions, the secrecy rate can be 
written as 

1 (1 + lar) 

Rs = i: log2 7— T (21) 

2 (l + 7ae+7re) 

where 

^ar ^ ^ ilar^a v^^j 

-iae = Patf Hf,(H5eT',Q,fe,Tf H^ + a^Y^YLaeta (23) 

7,e = p.hf,(H,eT;Q,,,Tf Hf, + d^Ij-^h^e, (24) 

and we aim to find the joint optimal solution for the jamming beamformers T'^, T'j^, the covariance 
matrices Q^;,', Qzas and the transmit power vector p = \Va:Vr^Va2-,'Pb^ in order to maximize the 
secrecy rate Rs. 



We will first consider optimizing the jamming beamformers and covariance matrices. For T'^ and Qzb'^ 
the problem of minimizing the SINR at Eve jae can be written as 

min tfHf,(H,eT',Q,,,TfHg + a2l)-iH,eta (25a) 

s.t. tr(Q,feO < pb, hfe,T; = 0. (25b) 



Although problem p5| ) can be formulated as a semidefinite program (SDP) that can be solved efficiently 
(see Appendix |A]), we can not directly obtain an analytical solution that is useful for optimizing the global 
power allocation. Therefore, we will make use of the following lemma, a proof of which is provided in 
Appendix IB] 



Lemma 1: The covariance matrix Q^;,' that minimizes ( |25a| ) is rank one. 

According to Lemma [1] we know that a one-dimensional jamming signal is optimal for the case of 
single data stream transmission: T'^ = t'j^. Under the constraint that h^rff^ = 0, and defining G^ as an 
orthonormal basis for A/'(h{,r), the jamming beamformer from Bob can be written as t'^ = G^Cf,, for 
some unit-length vector Cf,. Eq. ([23]l becomes 

7ae = Patf Hf,(pbHbet;tf H£ + a'^iy^Haeta 
- P^ (^H^H^ . _ tfHg,Hfeet;tfHgH,eta \ 

where the second equality holds due to the matrix inversion lemma |32|. The optimization problem is 



equivalent to maximizing the second term in (26 1, which can be formulated as 



m_ax Hfa'j . J^M-r^ ^_ s.t. Cj Cfe = 1 (27) 



^^ c^(f^I + BfB,)c, 
where a;, = G^^H^Hagta and B^ = H^eG^. The maximum value of the Rayleigh quotient in (27l is 
the largest generalized eigenvalue of the matrix pencil (af,a^, |^I + B^Bfe), and the vector that achieves 



it is the corresponding generalized eigenvector |33|. Since a^a^ is rank one, the solution can be written 



as 






'Pb' 



and 7ae becomes 



a" V " "^ " Pb' 



^ ( tf Hf,H,eta - af (— I + Bf Bfe)-iab) . (29) 
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Similarly for the second phase, the SINR for Eve is rewritten as 

\-,H-rT i-' i-'H-riH-L^ 



"'"'^•^''^^ tf(£l + H^,H..)t'J- ^''^ 



Using the same method as in (26l-(28l, Alice's jamming beamformer is given by 



" ^ 11(^1 + B^B„)-ia,|| ' ^'^ 

where G-^ is an orthonormal basis for 7\A(h^Hafc), a^ = Gf ^H^hre,Ba = HaeG-^, and jre becomes 

7.e = ^ f hf,h,, - af (^I + Bf B,)-ia.) . 

Next we find the power allocation that maximizes the secrecy rate. Note that the jamming beamformers 
are not independent of the jamming power, and thus we need to jointly optimize over both quantities. In 
general, ( [2T] ) is not convex with respect to p, so instead we maximize the following lower bound for Rs : 

1 'y 1 Ih t |2 
-■- Jar '- \'-'-ar^a\ 

-Rs(p) > - log2 -— r = - log2 2 / \ ^'^^^ 

2 (1 + Jae + Ire) 2 Cr^5r(p) 



^(P) =Pa^ +Pb ^ +Pa ^PrPa2 (33) 

2 

p-1 = tf Hf,H,eta - af (— I + B^B6)-iab (34) 

Pb 
2 

p-:^ = hf,h,e - af (-1 + Bf Bj-^a,. (35) 

Pa 

Over the range of practical transmit powers, ph and pa2 can be accurately approximated as linear functions 



of pi, and pa2, which we denote by pb = cipb + C2 and pa2 = csPa2 + C4. Note that according to (34i, 



as ph increases, the second term can only increase in size, which means p^^ decreases, and hence p^ 
increases, which implies that ci is positive. As pt approaches zero, the second term approaches zero, but 
the first term is non-negative, so that implies that C2 > 0. Thus ci and C2 are both positive constants. 
Similarly, we can see that C3 and C4 in ( |35] ) are also positive constants. 

Using this approximation, the rate maximization problem under a global power constraint P becomes 



one of minimizing g{p) in (32 1: 



mill Pa^ +Pb ^ +Pa ^PrPa2 (36a) 
P 

S.t. Pa + q Vb < -P + C]'^C2, Pr + C^^Pa2 < P + C^^C4 (36b) 

Pa|harta|^ =Pr|hrfeP (36c) 
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where (36b]) is derived from Pa + Pb ^ P and p,,. + Pa2 < P, and ( 36c| ) is the optimal power adjustment 



for the two hops used to guarantee that ^ar = Irh- The optimization problem stated above is in the 
standard form for Geometric Programming (GP) problems, with ( |36a[ ) and ( |36b[ ) as posynomial and 
( |36c[ ) as monomial constraints. GP problems are a class of non-linear optimization problems that can 
be readily turned into convex optimization problems, and hence a global optimum can be efficiently 
computed |34| . If individual power constraints are employed, we can also use GP to solve the following 
similar optimization problem: 

mill Pa^ +Ph'^ +Pa ^PrPai (37a) 

p 

S.t. Pa<Pa, Pr<Pr (37b) 

qVfe < A + q ^C2, C^^Pa2 < Pa + C^^Ci (37c) 

Pa Ihartal"^ = Pr\Kb\'^ ■ (37d) 



Remark 1: As discussed in the beginning of this section, we choose the principal generalized eigen- 
vector of the pencil (I + 2|h^har,I + %H^Hae) as the information signal transmit beamformer ta- 
However, the allocated power pa is unavailable before the optimization algorithm starts. Therefore, 
iterations will be needed for computing the beamformers, initialized with pa = P, where P is the 
maximum transmit power. Based on our numerical experiments, the algorithm usually converges with 
very few iterations, and introduces little complexity to the overall algorithm. 

Remark 2: For the case where H^e does not have full column rank, i.e., Na > Nf., an alternative 
would be to choose t^ to lie in the null space A/'(Hae). This beamformer will in general be different 
from the one we propose, and will result in a solution where Eve will not receive any information signal 
in the first phase and hence the jamming from Bob is not necessary. However, based on our numerical 
experiments, the solution we propose yields a larger secrecy rate. This is mainly due to the fact that 
although ta may allow a small amount of information leakage from Alice to Eve, the rate improvement 
in the information channel outweighs that of the wiretap channel, given the cooperative jamming support 
from Bob and the optimized power allocation. 

2) Minimum transmit power with fixed secrecy rate: The problem of minimizing the transmit power 
under a certain fixed secrecy rate is similar to the problems discussed above. We still choose jamming 
beamformers that lie in the subspace orthogonal to the intended channels. As before, for the first phase 
we will have t'^ = pbG^^eb, where ph is a scalar that maintains the unit norm of t'^. We aim to minimize 



12 

the norm of e;, under a fixed target secrecy rate Rq. According to (27i and (32 1, the problem can be 
formulated as 

where f{Ro,t'^) is a function of Rq and t^ independent of t'^,. The solution is again seen to be the 
generalized eigenvector of the pencil (a^a^, |-I + B^Bfc) corresponding to the largest eigenvalue. Since 

I (^I+BfBb)-ia6 

it is a rank-one Hermitian matrix, the result can be explicitly presented as t, = Grr---^ — ^ -, 

f ^ f b b ||(^I+BfBO-ia,||' 

which is the same result as in the rate maximization problem. Similarly, for the second phase, we also 



have the same beamformer as (31 1. Considering the transmit power of all the nodes, we can now formulate 



the optimization problem under the global transmit power constraint as 

IhartaP 

min mayi{pa+Pb,Pr+Pa2) s.t. g(p) < ^°^ "^ (39) 



where g{p) is given in (32 1. This is also a GP problem. To minimize individual transmit powers, (39 1 



should be rewritten as m.mpm.ax{pa,Pb,Pr,Pa2) instead. 

B. Multiple data stream relaying 

Since ECSI is known to the relay network, Alice and the Relay can utilize certain beamformers to 
perform multiple data-stream relay transmission and reduce information leakage to Eve as well. The 
GSVD has been employed for the traditional MIMO wiretap channel pjf , and it operates by dividing 
the channels from the transmitter to the intended receiver and the eavesdropper into a set of parallel 
subchannels. 

Definition 1 (GSVD Transform): Given two matrices Hi G C^-''^" and Ha G C^^'^^- and k = 
ranfc{[Hf ,Hf]^}, there exist unitary matrices U € C^-^^-, V G C^=^^<= and * e C^»^^°, and a 
non-singular upper-triangular matrix R E ^J'^^ such that 

U^Hi* = Si[R,OfcxJV„-fc], V^Hs* = S2[R,0fcxiv„-fc] 

where Si € M.^-^^^, S2 G M.^--^'^ are nonnegative diagonal matrices with Sf Si + 8^82 = Ifc, the diagonal 
elements of (8^81)2 are ordered as < si_i < • • • < si^^, and the diagonal elements of (81^82)2 are 
ordered as S2,i > • • • > S2,fc > 0. 

It has been shown ||7| that, for the standard Gaussian MIMO wiretap channel, using the GSVD-based 
beamformer 



IR- 



R 1 
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to transmit the desired signals along dimensions where si_j > S2,i achieves the secrecy capacity in the 
high SNR regime with uniform power allocation. In this section, two transmission strategies based on 
the GSVD will be investigated for the two-hop relay channel. In the first strategy, each transmission 
phase is treated as a standard wiretap channel, and Alice and the Relay will use GSVD-based transmit 
beamformers in the first and the second phase, respectively, without any cooperative jamming from 
inactive nodes. In the second strategy, a cooperative jamming scheme is proposed in which Bob and 
Alice also transmit jamming signals based on the GSVD transform in a reverse manner. 

1) Simple GSVD-based relaying: To begin, we consider the case where GSVD-based beamforming is 
used without jamming. According to Definition 1, the MIMO channels in phase 1 and phase 2 can be 
decomposed as 



H„ 



UaSar[R-a,0 



sxNa 



1* 



H 



Hae — VaSae[Ra, 



SXNa — S 



* 



H 



Urb — UrSrfe[Rr,0, 



sxN^~s 



* 



H 



XX|^g — V ^O^g XV-^. vl 



sxNr.—s 



* 



H 



where s = mm(ran/i;{[H^,H^]'^},ranA;{[H^,H^]^}), representing the maximum possible number 
of data streams. Alice and the Relay transmit signals with the following two beamformers respectively, 



R-a 



On„-sxs 



*r 



Rr 



-1| 



R^ 

On,.-sx 



(40) 



Proposition 1: When ( |40| ) is used for transmit beamforming, the secrecy rate under the simple GSVD- 
based relaying scheme can be expressed as: 



1 mm 

Rgsvd = ^ log2 



{rc=i 



1 + 



nti 1 + ^ 



Pa.iS] 



^I\i=i (1 + ^2['|R;;;''i[|2J| 



+ 



P-r.iS,. 



(41) 



where pa^i and pr^i are the transmit power for the ith parallel channel from Alice and the Relay, 
respectively. 

The proof of Proposition [T] is given in Appendix |C] Next, we will investigate the power allocation 



for the above transmission scheme. Maximizing the rate in (41) is generally a nonconvex optimization 



problem. However, applying the single condensation method for GP |35 1, the posynomial in the numerator 



of (41) can be accurately approximated as a monomial, and we can still solve this nonconvex problem 



through a series of GPs. 
Lemma 2: Let 



n/*(^--o=n 1+ 



Pa,iS„ 



i=l 



4 = 1 



n-2||R-M|2 



(42) 
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We have 

n/.(p.,.) > n/i(P„,.) = n (^r {-^§^^ 

i=l i=l i=l V"l.»/ V«2,*CT \\tta II ^ 

where ai^j,a2,i > 0. The inequaUty becomes an equaUty when ai,i,a2,i satisfy 



1 



Pa,i dfi{pa,i 



ai,i 



a2,; 



fi{Pa,i) ' "^'' /i(Pa,i) <9pa,j 

in which case 11^=1 fi{Pa,i) is the best local monomial approximation of ni=i fi{Pa,i) near pa,i 
Proof: We can rewrite ni=i fiiPa,i) as 

n/ -, , Pa,iSar.i \ TJ I ^ , Pa,i^ar,i \ 

^ + 9MT^-lll9 =11 "l.i + "2,i ^||„-i||^ 



(43) 



(44) 



i=l 



rT2||R-l||2 



> 



i=l 



n 



1 

ai < 



Pa.iSa 



a2.iO- 



IR"H|2 



(45) 



(46) 



where (46 1 holds according to the arithmetic-geometric mean inequality. Noting that ai^j and a2,i are 



both positive coefficients and ai^i + 02, j = l,Vi, the proof of equality is straightforward by inserting 



(|44|| back into llt=ifiiPa,i)- 

Similarly for the second phase, given the posynomial 



n5'(p-.*)=n^+ 



2 

PrA^rb'i 



i=l 



4 = 1 



we have the approximation 



n5.(p..)-n(^) 

1=1 j=i ^ ' ' 



ft,. 



where 



/3: 



l,i 



,/3; 



^2||R-l||2 I ' 
/ Pr,iSrb,i 

Pr,i dgi{pr,i) 



'2,i 



9i{Pr,i) "^^'' S'iCPr.i) 5pr-,i 

The approach corresponding to these results is outlined in the following algorithm: 
Algorithm 1: Single condensation method for power allocation 
Initialize p)^j and p;/, i = {1, • • • , s}. 
For iteration k: 

1) Evaluate posynomial fi{p\i ) and gi{p\.^ ), according to (42 1 and (47 1. 

2) Compute a^'^) and /S^*^); 



(47) 



(48) 



(49) 



a 



(k) 

l,i 



1 



(fc) 



1 



a. 



(k) 

2.i 






Ak-l) 



hij^r') diC 



(fc-1) 



9iipZ ^^) 
' P2,i - 



Jfe-1) 



P;,i ' 9gi{p 



(fc-i). 



(50) 



Jfc-i) 



9i{PrA ') dpi- 



(fc-1) 
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3) Condense posynomials fi and gi into monomials fi and gi, according to (43 1 and (48 1 

4) Solve the GP 



2 2 



minmax lY[Mpa,i) \llg^{Pr,^) ' m 1 + ,2||R-1||2 + ^2|iB-'l||2 ^^la) 

s s 



i=l i=l 



5) Apply the resulting p^ ■ and p^ ■ to step 1 and loop until convergence. 

The GP problems in this successive optimization method can be solved using interior-point methods 



with polynomial-time complexity |36|, and it has been proven in |35| that the solution obtained using 
successive approximations for the single condensation method will efficiently converge to a point satisfy- 
ing the KKT conditions of the original problem, and the global optimum can consequently be obtained. 



Note that (BTab is refered to as a generalized posynomial |[36i since it is formed from posynomials using 



a maximum operation, and can be easily converted to the standard posynomial form as 

^ + 2II-D-1II2 + 2IIP-1II2 *^^^^) 

s s 

S.t. ^Pa,i<P, ^Pr,i<P (52b) 

i=l 4=1 

s s 

llUpa,^r^ll-^<l, llUPr,i)~^l^^^<i- (52c) 

j=l j=l 

2) GSVD-based PCJ: A GSVD-based, partial cooperative jamming scheme is proposed in this subsec- 
tion. In this case, Alice and the Relay will still use the same transmit beamformers as in the case without 
cooperative jamming. Since Bob and Alice are normally inactive in phase 1 and phase 2 respectively, 
they can act as temporary helpers to help improve the secrecy rate. As before, however, the power used 
for jamming must come from the total power budget of P in each hop. A GSVD-based beamformer for 
the jamming signal is used by Bob in the first phase, due to the assumption that ECSI is available. The 
GSVD is implemented in a reverse fashion, since Bob in phase 1 considers Eve as the intended receiver 
of the jamming and wants to avoid leaking interference signals to the Relay. Similarly in phase 2, Alice 



treats Eve as the intended receiver. The signal model for this scheme is given in Section II-B 

Performing the GSVD for the channels from Bob to Eve and the Relay according to Definition 1, we 
have 

Hfee = Uf,Sf,e[R-fe, Ofc,xAr,-fcJ*6 Hbr = VfeSfer[Rfe, Ofc.xTV.-fcJ^b 

Hae = Ua'Sae[R-a',Ofc^xAr„-fc„]*a' ^ab = Va'Sab[Ra' , Ofc^x Ar„-fc„]*a' 
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where k^ = ran/i:{[H^, H^]^}. Bob and Alice use the following jamming beamformers to implement 
the reverse GSVD 



Tl 



*ft 



ir;; 



-1 



R 



ONr-kbXkb 



T' 



*^ 



IR 



a' 



^.' 







Na-kaXka 



and, unlike the simple GSVD-based relaying scheme, there will be jamming energy present in the signals 
received by the Relay and Bob. For Eve, the received signal is given by 



H„f,T„D 



ae -*- a^-'a 



X i^g -I- ^ 1^ f 



Z + 



HgZ + ne. 



H5eT'j,z; + riei 

Employing the above jamming beamformers, the mutual information between Alice and Bob is 

Id = mill I ^ log2 det(I + H„,T„Q,„Tf Hf^Qr^), ^ logs det(I + H,feT,Q,,Tf Hf^Qr^^) 
where E(z;zf ) = Q.^' ^(z^zf ) = Q,,,-, and 

Q^, = E[(Hbr%4 + n,)(Hb,T;z'b + n,)^] = Hb^T^Q^fe-Xf H^ + a^I 



The mutual information at Eve is 



1 



le = mill <j ^ log2 det(I + H„,T,Q,„Tf Hf^Q-^), - logs det(I + HeQ.Hf Q-^^) 



where 



Qf 



HbeT'Q,fe,TfHg+CT2l 







UaeKQ^a'K^ng + a^I 

To maximize the secrecy rate, we then have the following optimization problem: 



max 



R. 



PC.} 
'gsvd 



s.t. tr(Q^a + Q,b') <P, tr(Q,, + Q,a') < P 
where R^^i = Id — le and P is the global power constraint. 



(53) 



(54) 



(55) 



Remark 3: The secrecy rate in this case does not have a form similar to (41 1, and finding the optimal 
power allocation for this case is generally intractable. Therefore, we will use Newton's method initialized 
with the optimal point from the GSVD-based relaying algorithm. Though this may not find the global 



optimum, we can at least gain insight into this strategy. The global power constraints in (55 1 are set for 
fair comparison with the case without cooperative jamming. 
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IV. Secure relaying with unknown ECSI 

In this section, we assume that ECSI is unknown to the relay network. Thus, AUce and the Relay can 
no longer use beamforming methods like those based on the GSVD to selectively transmit information 
away from and jamming signals towards the eavesdropper. However, cooperative jamming can still be 
used to improve the secrecy of the information in the two-hop network. As described below, the approach 
we take to achieve this goal is to first meet a fixed target rate for the relay link, and then allocate all 
remaining resources to wide-area jamming, while guaranteeing that the jamming signal has no impact 
on the desired information. 

We propose a cooperative jamming strategy in which the signal space is divided into two orthogonal 
subspaces, an information subspace and a jamming subspace. Both PCJ and FCJ approaches can be 
applied in this scenario. For PCJ, any available jamming power will only be allocated to information 
transmitters, while Bob (phase 1) and Alice (phase 2) remain inactive. For FCJ, both the transmitter 
and the temporary helpers can perform cooperative jamming in the jamming subspace, which will allow 
the legitimate receivers to use beamforming to reject interference from this subspace. Note that when 
using FCJ, cooperative jamming requires the receiver to broadcast the jamming subspace so that the 
interference can be aligned at the desired receiver without a loss of information. Although Eve may also 
be aware of this subspace, she can not remove the jamming signal since she sees different channels from 
the transmitters and jammers. 

In phase 1, assume spanjHar} = span{r]i, r]2, . ■ ■ , rik,rik+i, • ■ ■ , Vm}, where k is no greater than the 
maximum possible number of data streams, and r]i,rj2 . . . ,r]m form an orthonormal basis. The information 
and jamming subspaces are defined to be Si and J'l, respectively, where Si = span{rji, r]2, ■ . ■ , %} and 
J'l = S^. Assuming the receive beamformer matrix at the Relay is W,. = [r/i, ry2, • • • ,?7fc], the signal 
received by the Relay is 

y^ = W^[Uar{TaZa + T^z'J + Ubr%4 + n^] = Uar^a + Hr- (56) 

where Har = W^HarTa, z^ is the information signal vector transmitted by Alice with covariance 
Qza, z'(j and z'^ are jamming signals transmitted by Alice and Bob, with covariance matrices Qz'a and 
Qz'b, respectively. The transmit beamformers are chosen such that HarTaZa G Si, and Ha^T'^z'^ G 
Ji,HferT'jjZ'^ G Ji. The signal received by Eve in phase 1 is 

yel = HaeiTaZa + T^z'J + Ube%z'^ + U^l = HaeZa + Hel (57) 

where Hei = HaeT'^z'^ + H^T^z'^ + riei. 
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In phase 2, signal ^2 and jamming J2 subspaces are chosen from spanjHrft}, and similar to phase 
1, the signals at Bob and Eve are 

yfe = Wf [H,6(T,Z, + T'X) + nabT'^2<2 + nfe] = HrfoZr + Ub (58) 

ye2 = HreiTrZr + T[,Z^) + JlaeT'^2^'^2 + ^62 = H^eZr + ne2 (59) 

where Urb = W^H^tTr, H^e = H^eTr and ne2 = H,,eT'j,z'^ + HaeT'^2Za2 + '^e2> Zj, is the information 
signal transmitted by the Relay with covariance Q^,.. z'^2 ^^'^ '^'r ^^ jamming signals transmitted by 
Alice and the Relay, with covariance matrices Qz'a2 and Qz'r, respectively. As before, the beamformers 
T, and T; force H^,bT,.z,. e ^2, and UrbT'X G J2, llabTa2^'a2 G ^2- 

The cooperative scheme outlined in this section involves the allocation of power and the number 
of dimensions for the information and jamming subspaces. If the MIMO channel is rich enough, more 
dimensions allocated to the signal subspace increases the amount of power available for jamming, but leads 
to a smaller dimensional jamming subspace for both transmitters and cooperative jammers. More antennas 
for Eve usually requires a higher dimensional jamming subspace, especially when ECSI is unknown to 
the transmitters. One of the advantages of FCJ in this case is that in addition to the pre-assigned jamming 
subspace of dimension Na — k (for phase 1), the helpers provide jamming support in additional dimensions 
due to the fact they have different channels to Eve. Taking the transmission in phase 1 as an example, 
assuming k dimensions are assigned to the information subspace, the jamming subspace seen from Eve 
will be greater than Na — k. In particular, Na — k< dim{span{'iiaeT'a} H span{'H.be'^'b}) — ^(-^a — k). 

Therefore, the tradeoff between power and allocation of the jamming subspace dimension needs to 



be considered. In this case, we propose to use an approach similar to that in |24| and minimize the 
product of the power allocated to the information signal and the dimension of the information subspace, 
{Pa +Pr)k, such that the fixed target rate for the relay transmission is achieved. We then allocate all the 
remaining dimensions and power for jamming. Since the ECSI is not known, the jamming power will 
be uniformly distributed among all available dimensions at the transmitters and cooperative jammers. 
Assuming the target rate for the relay transmission is Rt, we have the following FCJ algorithm: 
Algorithm 2: FCJ with unknown ECSI 

1) Initialize Svd(Har) = Var'^ar^ar ^nd Svd(H.rb) = Vrb'^rb'^rb- 

2) While i < s 

. Let W^ = Uar[:, 1 : i], Wb = U,b[:, 1 : i], T^ = \ar[-, 1 : i], T, = V,fe[:, 1 : i]. 
. Let T'„ = Yar[-;i + 1 : Na], % = Yrb[:,i + 1 : Nr]. 
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Let SvdiW^Hbr) = Vbr^br^jf^, % = Vbr[:,i + 1 : N^], and Svd{Wj^Uab) = Vab^ab^^b' 

Solve the following problem 



pi^) = min tr{Q,a), pI'^ = min tr(Q^r) 

s.t. - log2 det(I + ^Har-Q^aHf^) = Rt, 77 loga det(I + - 
2 a'^ Z a' 



where the water filling algorithm is used to determine Q^^ and Q^r- 

3) Find k = argmirij [pi + Pr \ ■ i, and determine all beamformers for the resulting k. 

4) Allocate pi to diag{Qza}, and pi- to diag{Qzr} using water filling. 

5) Uniformly allocate P — pi' to diag{Qz'a, Qz'b}, and P — pi- io diag{Qz'r, Q2'a2}- 

The PC J algorithm in the unknown-ECSI case is similar to that for FCJ, except that jamming support 
will not be provided by Bob (in phase 1) and Alice (in phase 2), and thus the beamformers T'^ and T'^2 ^^ 
step 2 will not be used. In step 5, when the necessary amount of power for information signals is assigned, 
all remaining jamming power will be used by Bob and Alice in phase 1 and phase 2, respectively; i.e., the 
power P — pa and P — pi- will instead be assigned to diag{Clz'a} and diag{Qz'r}- In either approach, 
the optimization problem in step 2 can be solved with a simple line search. If the minimum rate Rt 
can not be achieved with the available power, the link is assumed to be in outage. In this algorithm, we 
assume that the Relay uses the same information dimension as Alice, as discussed in Section |ll] However, 
using different information dimensions for the two phases with a more complicated coding scheme may 
also be an interesting case to consider for future work. 

V. Numerical Results 

In the following simulations, the elements of all the channel matrices are assumed to be i.i.d. complex 
Gaussian. As shown in Fig. [2J Alice, Bob, the Relay and Eve are assumed to be located at (—0.5,0), 
(0.5,0), (0,0), {dx, —0.5) respectively, where distances are expressed in kilometers. We adopt a simple 
transmission model in which the standard deviation of each channel coefficient is inversely proportional to 
the distance between two nodes. We assume a path-loss coefficient of 3, and the same background noise 
power (T^ = — 60dBm at all nodes. All results are calculated based on an average of 3000 independent 
trials. 



For the known-ECSI case discussed in Section III we examine the performance of the following three 



schemes: PCJ for single data stream relaying (Section III-A| ), simple GSVD-relaying (Section III-Bli 



and also GSVD-PCJ (Section III-B2| ) for multiple data streams. For the unknown ECSI case discussed in 
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(-0.5,0) (0,0) (0.5,0) 

I km 



(dx,-0.5) ' E 



Fig. 2. Simulation scenario showing locations for Alice(A), Bob(B), Relay(R) and Eve(E). 



Section IV both the FCJ and PCJ approaches are simulated. For each of these schemes, we also examine 
the impact of both global and individual power constraints. For the case of individual power constraints, 
we assume the total transmit power to be evenly distributed to all transmit nodes. Furthermore, in order to 
examine the performance gain of the proposed cooperative jamming schemes and optimization algorithms, 
we also investigate cases using uniform power allocation, as well as cases involving conventional relay 
transmissions without jamming. 

The secrecy rate as a function of transmit power is shown in Fig. [3] for a case with known ECSI, 
where Alice and Bob both have four antennas, and the Relay and Eve each has one. Eve is assumed to 
be located closer to the Relay at (0, —0.5), which (as will be seen in the next example) is usually the 
worst-case assumption for the relay link. This will be the default assumption unless otherwise specified. 
Compared to traditional DF relaying, the PCJ schemes provide a significant improvement in terms of 
secrecy rate in the medium and high SNR regime. The benefit of having the flexibility associated with a 
global power constraint over fixed individual power constraints is clearly evident. Also, the performance 
gain of using geometric programming for power allocation is obvious, compared to the uniform power 
allocation scheme. We can also see that even the conventional relaying scheme is better than PCJ schemes 
with individual or uniform power allocation when the transmit power is low. This is because, with a less 
flexible power adjustment, a fraction of power that could have brought higher secrecy rate if used for 
data transmission is wasted on jamming signals. This illustrates the importance of an efficient power 
allocation if cooperative jamming support is applied. 

Fig. |4] presents the impact of Eve's location on the transmit power fraction for both the information 
and jamming signals, assuming that the global transmit power is limited to lOdBm. Unlike the settings in 
Fig. [3} Eve has four antennas in this scenario, which provides her with increased eavesdropping abilities. 
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Fig. 3. Secrecy rate vs. transmit power for Na — 4,Ni, — 4, N^ — 1, Nr = 1, and Eve located at (0, —0.5). 



In this case, we plot the secrecy rate performance as Eve moves from (—1,-0.5) to (1,-0.5). The 
secrecy rate is smallest when Eve is at the midpoint (0, —0.5), and increases in either direction away 
from (0, —0.5). Note also that the fraction of the transmit power devoted to jamming also decreases as 
Eve moves away from the midpoint. This behavior is due to the fact that, when Eve is closer to either 
Alice or Bob, most of her information about the desired signal comes from only one of the hops, due 
to the fact that the other hop is farther away and can be effectively jammed with minimal power by the 
transmitter she is closest to. This is the primary benefit of the cooperative jamming support provided by 
PCJ. 

The performance of GSVD-based relaying without cooperative jamming and GSVD-based PCJ strate- 
gies, where the relay link transmits multiple data streams, is shown in Fig. [5] Here we see that cooperative 
jamming with global power allocation provides considerable gain in secrecy rate over other schemes. 
However, the use of individual power constraints significantly degrades the benefit of the jamming signals, 
although it still approaches and even surpasses the performance of GSVD-relaying with optimal power 
allocation when the transmit power is higher. In addition, we also see the benefit of Algorithm 1 for 
power allocation in the GSVD-relaying scheme, compared with using simple uniform power allocations. 

Finally we consider examples for the case where ECSI is not available. The transmit power P is set 
to be 15dBm in these examples. In Fig. [6l all nodes are equipped with four antennas, and the secrecy 
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Fig. 4. Secrecy rate and transmit power fraction vs. eavesdropper location for Na = 4,Nb — 4, A^'e = 4, A^,- — 1, global 
power constraint P = lOdBm, and Eve's location varies from (—1, —0.5) to (1, —0.5). 
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-A — GSVD-PCJ global power allocation 
-e — GSVD-PCJ individual power allocation 
-B — GSVD-relaying optimal power allocation 
H - GSVD-relaying uniform power allocation 
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Fig. 5. Secrecy rate vs. transmit power for Na — 4,Nt — 4, A'^^ ~ 4, A''^ = 4, and Eve located at (0, —0.5). 



performance is given as a function of the rate constraint at Bob. For purpose of comparison, a naive PCJ 



scheme that uses the criterion min(pa +Pr) (instead of min(pa +Pr)k as discussed in Section IV I is also 
simulated. It can be seen that if no jamming signals are used, there is little difference between the mutual 
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information at Bob and that at Eve, and thus we expect the secrecy of the message to be low. Similar 
to Fig. [3} the individual power constraint will reduce the secrecy performance due to the inefficiency of 
the power assignment. We can see that FCJ achieves a big performance gain compared with PCJ as Rt 
increases, since FCJ leads to a higher dimensional jamming subspace than PCJ, although they transmit 
with the same jamming power. In addition, the performance of PCJ begins to level off and even drop for 
high Rt, since more power is allocated to information signals, and the protection from eavesdropping is 
reduced. 

Fig. [7] provides a detailed look at how the number of eavesdropper antennas affects the performance 
of the different cooperative jamming schemes. In this case, we fix the target rate for relay transmission 
to be Rt = Ibps/Hz. Alice, Bob and the Relay are equipped with four antennas, and the number of Eve's 
antennas increases from one to eight. It can be seen that when Eve has only one antenna, little advantage 
is observed for FCJ since Eve only receives single-dimensional signals. However, as the capability of the 
eavesdropper increases (i.e. when Eve has more antennas), the relative gain of FCJ over PCJ increases, 
although the performance of all methods decreases. 



- FCJ global power allocation 
A - FCJ individual power allocation 

-e— PCJ 

e - PCJ naive 
— s — w/o jamming 




Fig. 6. Secrecy performance vs. rate constraint for relay link when ECSI is unknown, Na — 4:,Nb — 4, A^e 
Eve located at (0, -0.5), P = 15dBm. 
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Fig. 7. Secrecy performance vs. eavesdropper antenna number when ECSI is unknown, A^^ = 4, A^j, = 4, A^^ = 4, Eve located 
at (0, -0.5), fixed target rate Rt = Ibps/Hz, P = 15dBm. 



VI. Conclusions 

In this paper, we have proposed partial cooperative jamming (PCJ) and full cooperative jamming 
(FCJ) strategies for two-hop DF relay systems in the presence of an eavesdropper that can wiretap 
both transmission phases. Both single and multiple data stream transmission scenarios were considered. 
For single data stream relaying, the system design was conducted from the perspective of secrecy rate 
maximization and transmit power minimization. By adopting the zero-forcing constraint that the jamming 
signals be nulled out at the intended receivers, we obtained closed-form expressions for the jamming 
beamformers and the corresponding power allocation. For the case of multiple data stream transmission, 
we proposed a GSVD-based relaying scheme without jamming, as well as a GSVD-based PCJ scheme. 
The latter shows a significant performance improvement even though only a potentially suboptimal power 
allocation scheme is used. We also studied the secure relaying problem when the eavesdropper's CSI is 
unknown. Instead of maximizing the secrecy rate, a more reasonable relaying scheme with both PCJ and 
FCJ is proposed in which a target QoS for the relay network is met, and only the remaining resources 
are used for jamming. These schemes are shown to provide large gains in terms of the difference in the 
mutual information between the desired receiver and the eavesdropper. In particular, FCJ is shown to be 
a better choice when the eavesdropper's CSI is unavailable since the ability to exploit additional jamming 
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subspace dimensions is preferable when the transmitters possess no information about the eavesdropper. 



Appendix A 



SDP FOR PROBLEM (25 ) 



Let Qzb' = T^'bQzb'T^'b^ ^ where T'^ is normaUzed such that ti{Qzb') = tr(Qzfe')- Problem (25 1 is 



equivalent to 



mm /i 

Qzb'bO,At>o 



s.t. t^Hf,(HfeeQ.fe'H^ + a^iy^Haeta < /" 



tr(Q 



zb' 



< 



Pb- 



Using the Schur complement p2| , constraint ( |60b| ) can be written as 

rH4.H 



^ 



■CT-H -t-J 



Haeta HfegQzb'H, + (T I 



^0. 



(60a) 

(60b) 
(60c) 



(61) 



Combining ( [61] ) with the trace constraint and the semidefinite constraints on Q^;,' , the equivalent problem 
becomes 



mm jj. 

s.t. tr{Q^b,) < Pb, Q^b' ^ 0, ^ > 0, hbrQzb' 



fJ- 






^0. 



(62a) 
(62b) 

(62c) 



Haeta HbeQzb'H-be + <^ ^ 

This is an SDP that consists of a linear objective function, a linear equality constraint, and a set of 



linear matrix inequalities (LMIs) |34|, and thus can be solved efficiently, and T[, can be obtained via the 
eigenvalue decomposition of Q^^' . 

Appendix B 
Proof of Lemma[T] 



Given any jamming beamformer T'^, ^25\ becomes 



mill fiQzb') s.t. tr(Q2t.) <Pb 

Qzb'bO 



where fiQ^b') = t^Hf,(H6eT'^Q.fe'Tf Hg + aH)-^IIaeta, and the Lagrangian of ^ is 
L(Q,b', A, *) = fiQzb') + A(tr(Q,bO - pb) - tr(*Q,t,0 



(63) 



(64) 
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where A > is the Lagrange multipUer associated with the inequaUty constraint ti{Qzb') < Pb, and 
$ ^ is the Lagrange multipUer associated with the inequaUty constraint Q^;,- >z 0. Next, we can obtain 
the necessary conditions for the optimal Q^b- by using the Karush-Kuhn-Tucker (KKT) conditions: 



ti{Qzb')<Pb, Qzb to, A>0, *^0 
tr(*Q,;,0 = 
X{ti{Q,b') - Pb) = 
$ - = AI 
where is obtained by differentiating f{Qzb') with respect to Q^;,', and is given by 

= -HfteT'b (HheT'fjQzb^T'fj Hf,g + cr I)" 'H.aetata'H.ae(^beT^bQzb'% Hj,g + CJ I)~ T'^Hfee- 

Since ta is a vector, it is obvious that is a rank-one negative semidefinite matrix. 



(65) 
(66) 
(67) 
(68) 



For the case that A = 0, according to ( [68] ), we have = $. Since has a negative eigenvalue, $ 
will also have a negative eigenvalue, which contradicts the fact that $ is positive semidefinite. Thus A 



can only be positive. For A > 0, according to ([68]), we know that $ — is a positive definite matrix. 
Therefore, $ has at least A^ — 1 positive eigenvalues, i.e. rank(^) > A^ — 1, in order to keep ^ — @ >~ 0. 
Assuming Ai($) and Xi{Qzb'),i = {1> 2, . . . , A^} are eigenvalues of ^ and Q^;,', respectively, in non- 
increasing order, and due to the fact that $ and Q^f are both positive semidefinite matrices, we have 



yN 



tr{^Qzb') > J2i=i '^ii^)^N-i+iiQzb')- Combining this observation with ([66]), we also have 



N 



^Xi{^)XN^i+i{Qzb') = 0. 



(69) 



i=l 



Thus we can conclude that rank{^) ^ N, since otherwise all eigenvalues of Q^b' are zero and no 
jamming signals are transmitted. Combining this conclusion and the observation that rank{^) > N — 1, 
we can conclude that rank{^) = N — 1. Therefore, according to ( [69] ), we have Xi{Qzb') > and 
Xi:^i{Qzb') = 0, which indicates that rank{Qzb') = 1, and the proof is completed. 

Appendix C 
Proof of Proposition [1] 



According to the signal model given in Section II-A the signals received by Eve during both phases 
can be combined together as 



-H-ae ^a^a 



3. X7"g -I- T> J— ./t" 



z + 



nei 

ne2 



HeZ + He 



(70) 
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where D^ = diag{y/pk^} and IE(zz-'^) = Q^ = I. 

Using the transmit beamformers in ( |40| ), and denoting E(zaZ^) = Q^^ = diag{pa,i,- ■ ■ ,Pa,s}, 
E{zrZ^) = Qzr = diag{pr^i, ■ ■ ■ ,Pr,s}, the mutual information between AUce and Bob is 



- log2 det(I + -^UarTaQzaT^^Har), ^ log2 det(I + ^ 



Id = min <! ^ log2 det(I + ^H,,T,Q,,Tf Hf,), - logs det(I + ^H,bT,Q,,Tf Hfj \ (71) 



where 



^ log2 det(I + ^Ha,TaQ,aTf Hf,) = ^ logs det(I + ^S^.Q^aSf,) 
i log2 det(I + iH,i,T,Q„T?HS) = ^ logj det(I + J^S,,iq„Sft) 



For Eve, we have 



- logs det(I + -^UarQzallar), o ^^^2 det(I + ^ 



/e = min <j ^ logs det(I + -^Uar-QzaHar), o log2 det(I + -^HeQ.Hf ) \ (74) 



where 



^ logs det(I + -^HeQ.Hf ) = J logs det(I + ^Hf H^ 



^s - -- e / 2 °" ' cj' 
-log2det(I + - 

2 (T 



logs det(I + ^(Df Sf^SaeDa + Df Sf,S,eD,)) 



1 , "TT / - . Pa^i^ae.i Pr.i^re 



and according to the same secrecy constraints in (16 1, the secrecy rate (41 1 can be obtained. 



(72) 



1=1 \ M ' II / 



^"g^I 1 + ^2||R-1||2 + ^2||R-1||2 ' (75) 
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